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LONG-TERM X-RAY SPECTRAL VARIABILITY IN SEYFERT 1 GALAXIES 
A. Markowitz\ R. Edelson\ S. Vaughan^ 

ABSTRACT 

Direct time-resolved spectral fitting has been performed on continuous RXTE monitoring of seven 
Seyfert 1 galaxies in order to study their broadband spectral variability and Fe variability charac- 
teristics on time scales of days to years. Variability in the Fe Kq, line is not detected in some objects 
but is present in others, e.g., in NGC 3516, NGC 4151 and NGC 5548 there are systematic decreases 
in line flux by factors of ~2-5 over 3-4 years. The Fe Kq, line varies less strongly than the broadband 
continuum, but, like the continuum, exhibits stronger variability towards longer time scales. Relatively 
less model-dependent broadband fractional variability amplitude (Fyar) spectra also show weaker line 
variability compared to the continuum variability. Comparable systematic long-term decreases in the 
line and continuum are present in NGC 5548. Overall, however, there is no evidence for correlated vari- 
ability between the line and continuum, severely challenging models in which the line tracks continuum 
variations modified only by a light-travel time delay. Local effects such as the formation of an ionized 
skin at the site of line emission may be relevant. The spectral fitting and Fyar spectra both support 
spectral softening as continuum flux increases. 

Subject headings: galaxies: active — galaxies: Seyfert — X-rays: galaxies 



1. INTRODUCTION 

There is substantial evidence to indicate that the cen- 
tral regions of Seyfert 1 Active Galactic Nuclei (AGN) 
consist of a geometrically-thin, optically-thick accretion 

■ disk surrounding a supermassive (10^~^Mo) black hole. 
. In the leading models, the inner disk is surrounded by a 
' hot, optically-thin corona that Compton-upscatters soft 

photons thermally emitted from the disk to produce an 

i; X-ray power-law continuum over the range 1-100 keV, 

' though the geometry of the corona remains uncertain (e.g., 

[ Haardt, Maraschi & Ghisellini 1994). Furthermore, some 

' optically thick material, such as the accretion disk, inter- 

■ cepts and reprocesses some of the X-rays, producing the 
, so-called Compton reflection hump above '--'10 keV and a 
• strong Fe Kq fluorescent emission line at 6.4 keV (Guil- 
; bert & Rees 1988, George & Fabian 1991). Though the 
. exact configuration of accreting material throughout the 
' central engine is not known, the Fe Kq, line, as a primary 

signature of the reprocessing, serves as a potential tool to 
trace the distribution of accreting gas. 

Single-epoch spectral fitting has revealed a variety of 
profiles for the Fe Kq line. In some Seyfert 1 galaxies, the 
line is seen to be relativistically broadened (e.g., FWZI up 
to 100,000 km s~^ in the well-studied broad line profile of 
MCG-6-30-15, Tanaka et al. 1995) and redshifted. This 
can be explained in terms of a line originating in the in- 
nermost regions of the accretion disk, in close proximity to 
the black hole, and affected by strong Doppler and gravi- 
tational redshifts (Fabian et al. 1995). However, the ma- 
jority of Seyfert 1 galaxies observed with XMM-Newton 
(Reeves 2002) and Chandra (Yaqoob, George & Turner 
2001) feature a narrow Fe Kq component, modelled to be 
only mildly broadened, with FWHMs of <10,000 km s~^ 
(e.g., Kaspi et al. 2002, Nandra 2001, O'Brien et al. 2001, 
Yaqoob et al. 2001). For these objects, a line origin in 



the outer accretion disk, putative molecular torus, broad 
emission line region or narrow emission line region may be 
more relevant. 

Studying the time-resolved variability of the line and 
continuum components provides a complementary analy- 
sis to single-epoch spectral fitting. For instance, Fe Kq 
line emission originating in close proximity to the central 
black hole is expected to respond rapidly to changes in the 
illuminating continuum flux (Stella 1990, Matt & Perola 
1992). The simplest models predict that if the reprocess- 
ing time scale is negligible, then variations in the Fe Kq 
line should track continuum variations, modified only by 
a time delay equal to the light-crossing time across the in- 
ner disk. However, recent studies have not been able to 
support this picture. Investigations into the behavior of 
the broad Fe Kq line in MCG-6-30-15 by Iwasawa et al. 
(1996, 1999), Lee et al. (1999), Reynolds (2000), Vaughan 
& Edelson (2001, hereafter VEOl), and Fabian et al. (2002) 
have confirmed the presence of rapid Fe Kq line variability 
in this object, but have found the line to vary considerably 
less strongly than the continuum, and have found no cor- 
relation between continuum and line variations. Yaqoob 
et al. (1996) found the red wing of the Fe Kq profile in 
NGC 7314 to respond to rapid continuum variations, while 
the core did not respond as rapidly. Conversely, Nandra 
et al. (1999) found that the narrow core of the line in 
NGC 3516 tracked the continuum flux, while the broad 
wings displayed uncorrelated variability. Petrucci et al. 
(2002) used XMM-Newton to obtain two observations of 
Mkn 841 separated by 15 hours; while the continuum flux 
and shape stayed nearly constant between the two obser- 
vations, the narrow line virtually disappeared. 

This paper aims to study the broadband spectral vari- 
ability of Seyfert 1 galaxies on time scales of days to years. 
This includes quantifying the variability of the Fe Kq line 
and determining the relation between the line and the con- 
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tinuum flux. To this end, samples have been assembled 
that probe spectral variability on time scales of days to 
weeks (hereafter referred to as "short" time scales) and 
time scales of months to years (hereafter referred to as 
"long" time scales), using archival monitoring data from 
the Rossi X-ray Timing Explorer (RXTE). RXTE piovides 
the largest effective area at 6 keV of any current X-ray in- 
strument flown to date and thus offers the highest signal- 
to-noise of the Fe K region of the X-ray spectrum. The 
data reduction and sampling are described in §2. In order 
to track the line behavior, direct time-resolved spectral 
fitting has been systematically performed and fractional 
variability amplitude spectra have been constructed, pre- 
sented in §3. The observational results are summarized in 
§4. The physical implications are discussed in §5, followed 
by a brief summary of the main results in §6. 

2. DATA COLLECTION AND REDUCTION 

The goal of this project was to obtain uniform samples 
to survey line and continuum variability behavior. Con- 
tinuous RXTE archival monitoring data turning public be- 
fore 2002 December, when these analyses were performed, 
were used, as well as the authors' proprietary observations 
of three Seyfert 1 galaxies observed during Cycle 7. 

All RXTE data considered in this paper were obtained 
with the proportional counter array (PCA). 3 18 keV 
count rates were estimated with the online W3PIMMS 
tool, using mean 2-10 keV archival count rates and pho- 
ton indices obtained either from previously published fits 
(e.g., Kaspi et al. 2001) or the Tartarus database. The goal 
of time-resolved spectral fitting was to obtain at least ten 
high-quality spectra with at least ~50,000 photons in each 
time bin; sources that failed to meet this criterion were 
excluded. This resulted in eight observations consisting of 
continuous monitoring on short time scales (e.g., sources 
observed once every 1-2 orbits for ~1 ks over a span of 
~30 days) and an additional seven observations covering 
long time scales (e.g., sources observed once every 3 or 4 
days for several years). There were six objects common 
to both lists. This included two MCG-6-30-15 observa- 
tions on short time scales, an 8-day intensive scan hence- 
forth denoted as short-term MCG-6-30-15 (I) and a 53- 
day observation denoted as short-term MCG-6-30-15 (II). 
This also included two MCG-6-30-15 observations on long 
time scales, one spanning 326 days denoted as long-term 
MCG-6-30-15 (I) and a 2093-day observation denoted as 
long-term MCG 6-30-15 (II). Additionally, NGC 7469 was 
monitored only on the short time scale. Table 1 lists the 
date ranges and source parameter information. 

The PCA consists of five identical collimated propor- 
tional counter units (PCUs; Swank 1998). For simplicity, 
data were collected only from those PCUs that did not suf- 
fer from repeated breakdown during on-source time (PCUs 
0, 1, and 2 prior to 1998 December 23; PCUs and 2 from 
1998 December 23 until 2000 May 12; PCU 2 only after 
2000 May 12). Count rates quoted in this paper are nor- 
maUzed to 1 PCU. Only PCA STANDARD-2 data were 
considered. The data were reduced using standard extrac- 
tion methods and FTOOLS v5.2 software. Data were re- 
jected if they were gathered less than 10° from the Earth's 
limb, if they were obtained within 30 min after the satel- 

^ See http://lheawww.gsfc.nasa.gov/users/keith/rossi2000/energy_re 



nte's passage through the SAA, if ELECTRONO > 0.1 
(ELECTRON2 after 2000 May 12), or ii the satellite's 
pointing offset was greater than 0?02. 

As the PCA has no simultaneous background monitor- 
ing capability, background data were estimated by using 
PCABACKEST v2.1e to generate model files based on the 
particle-induced background, SAA activity, and the dif- 
fuse X-ray background. This background subtraction is 
the dominant source of systematic uncertainty (e.g., in 
total broadband count rate) in RXTE AGN monitoring 
data (e.g., Edelson & Nandra 1999). Counts were ex- 
tracted only from the topmost PCU layer to maximize 
the signal-to-noise ratio. All of the targets were faint 
(< 40 ct s~^ PCU~^), so the most recent versions of the 
'L7-240' background models were used. Source and back- 
ground spectra were derived over the time intervals listed 
in Table 1. These time bins were chosen to ensure ade- 
quate signal-to-noise in the Fe line in the time-resolved 
fits, but also to ensure adequate temporal resolution in de- 
terming the variability characteristics of both the contin- 
uum and the line. Response matrices were generated for 
each gain epoch using PCARSP V.8.0. 

3. ANALYSIS 

The time-averaged fits, described in §3.1, were per- 
formed with the goal of minimizing the number of free 
parameters needed to produce an acceptable fit. These 
model fits were then applied to the time-resolved data, as 
described in §3.2. The errors on parameters derived from 
spectral fitting are discussed in §3.3. The time-resolved 
spectral fitting is discussed in §3.4. Continuum and line 
variability amplitudes are presented in §3.5. Temporal 
cross-correlation analysis is described in §3.6. Finally, 
the fractional variability amplitude {Fyar) spectra are de- 
scribed in §3.7. 

3.1. Time-Averaged Spectral Fitting 

For all analyses, data below 3.5 keV were discarded 
in order to disregard PCA calibration uncertainties be- 
low this energy and also to reduce the effects of warm 
absorbers in sources (e.g., MCG-6-30-15, Reynolds 1997; 
NGC 3783, Kaspi et al. 2002; NGC 5548, Mcastro et al. 
2000). Data above 18 keV were discarded as the source 
count rate became too low at higher energies to obtain a 
good detection in each time bin. All spectral fitting was 
done with XSPEC v.11.0.2 (Arnaud 1996). Most long- 
term observations spanned multiple gain epochs, requiring 
simultaneous fitting of data sets. 

For all observations, the simple model of a power- law 
modified by cold absorption resulted in large residuals near 
6 keV due to the Fe Kq, line; the addition of a emission 
component proved significant at greater than 99.99% con- 
fidence in an F-test. Given the limit of PCA calibration 
(data/model residuals <2-3%'^, and the low energy reso- 
lution of the PCA (~1 keV at 6 keV), the PCA is largely 
insensitive to the detailed shape of the Fe Kq line. Both 
a single Gaussian and a Laor model component (diskline 
model for a maximally rotating Kerr black hole; Laor 1991) 
were tested to parameterize the line emission. More com- 
plex fits such as double Gaussians or a Laor component 
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plus a narrow Gaussian were either not required or did 
not yield sensible results; the single-Gaussian and Laor 
fits alone gave sensible results. For the Laor model com- 
ponent, the inclination angle was kept fixed at 30°, and 
the outer radius kept fixed at 400 gravitational radii; in 
no case did thawing either or both of these parameters 
measurably improve the fits. The emissivity index (3 was 
constrained to fall in the range 0.0-9.9. Laor models gave 
significantly better fits (at >99% confidence in an F-test) 
than the singk^-Gaussian parameterization in seven of the 
fifteen time-averaged spectra. An implicit assumption in 
these fits is that the emission regions in each source do not 
evolve substantially (i.e., major changes in the geometry 
of the emitting regions) throughout the duration of the 
observation. 

For all sources except NGC 4151, neutral absorption 
was kept fixed at the Galactic value. For NGC 4151, mod- 
elling the column density with a value much higher than 
Galactic resulted in a significant (>99.99% confidence in 
an F-test) improvement in the model fit. The column den- 
sity values obtained, given the PCA calibration, were in 
reasonable agreement with each other and with the value 
of 4.6±0.2 X lO^^cm-i in the spectral fit by Schurch & 
Warwick (2002). 

For all model fits except the short time scale observa- 
tion of 3C 390.3 and the Laor model fit to the short time 
scale observation of NGC 3516, there were residuals above 
10 kcV that signalled the need to incorporate a Compton 
reflection hump. Including a PEXRAV component to model 
reflection of an underlying power law component off opti- 
cally thick and neutral material (Magdziarz & Zdziarski 
1995) provided the necessary improvement to the fits, sig- 
nificant at greater than 99% c;onfidence compared to mod- 
els with just a power-law plus Fe K„ line, modified by cold 
absorption. In the PEXRAV component, the inclination an- 
gle of the reflector was kept fixed at 30° , and the power-law 
cutoff energy was kept fixed at 100 keV. All element abun- 
dances were kept fixed at solar; as noted by VEOl a change 
in Fe abundance caused the absolute values of spectral fit 
parameters to change slightly, but the variability pattern, 
the focus of this paper, remained unchanged. 

The basic model thus used, a power-law component plus 
PEXRAV (where applicable) plus a parameterization for the 
Fe K„ line as described above, all modified by neutral 
absorption, gave reasonable fits to the data. The excep- 
tional signal-to-noise in these timc^-averaged spcictra in- 
variably meant that systematic errors associated with the 
background modelling and spectral response of the PCA 
dominated the uncertainties. Consequently, finding a sta- 
tistically acceptable time-averaged fit (with xt ^ ^) ^'<^^ 
unlikely even if the model was representative of the in- 
strinsic spectrum of the source. The data/model residuals 
seen were consistent with the limit of the PCA calibration. 
Tables 2 3 summarize the time-averaged fitting results for 
both the single-Gaussian and Laor models, respectively, 
on both time scales. 

While higher-resolution X-ray instruments can detect 
one or both of the Fe Kq- broad and narrow components, 
the PCA is sensitive only to the combined flux of the two 
components. Recent XMM-Newton model fits can provide 
insight as to which component, the broad or narrow, dom- 
inates the RXTE view in each object. {Chandra-HETG 



can additionally constrain the narrow component but is 
mostly insensitive to the broad component.) For instance, 
the equivalent width of the broad component has been 
measured to be much greater than that of the narrow 
component in both MCG 6-30-15 (Wilms et al. 2001) and 
NGC 3516 (Turner et al. 2002). The present RXTE fits 
also indicated that the Fe lines in these two objects are 
among the broadest of the sample, consistent with the bulk 
of the line emission originating in the inner accretion disk. 
In constrast, the narrow component seems to dominate in 
NGC 4151 (Schurch et al. 2002), NGC 5548 (Pounds et 
al. 2003, Yaqoob et al. 2001) and NGC 7469 (Blustin et 
al. 2003), with a small upper limit or no detection of the 
broad component. The present RXTE fits for NGC 4151, 
for instance, indicated the narrowest line of the sample, 
consistent with emission far from the central black hole. 
The equivalent widths of the broad and narrow compo- 
nents in NGC 3783 are approximately equal (Blustin et 
al. 2002), indicating that emission from both the inner ac- 
cretion disk and a more distant location is relevant in this 
source. 

3.2. Time-Resolved Spectral Fit Components 

The above models were then applied to \\iq time-resolved 
data. Spectral fits to time bins spanning gain epoch 
changes required simultaneous fitting of the two data sets. 
The Fe K^, line was always detected in each time bin, as 
the addition of a line emission component was significant 
at >99% confidence in an F-test for nearly every time bin 
(in the single exception, the last time bin of the long time 
scale observation of NGC 5548, it was significant only at 
98.1% confidence). While fitting the time-averaged models 
to their respective time-resolved spectra, it was not practi- 
cal to thaw parameters siich as the Fe Kq line profile shape 
(line width in the single-Gaussian models; inner radius and 
emissivity index (3 in the Laor models) or line peak energy, 
given the PCA resolution and sensitivity. The F-test was 
used to determine which free parameters to include in the 
fit. The time-resolved spectral fits were repeated with each 
of photon index F, reflection fraction R (ratio of the re- 
flection normalization at 1 keV to the normalization of the 
power law) and line flux normalization V Ka frozen at their 
time-averaged values, and the values of total were com- 
pared via an F-test to determine which parameters could 
justifiably be allowed to be free and which should be kept 
fixed in the fitting, as listed in Table 4. There was no for- 
mal significant evidence (>99.7% confidence) of the need 
to thaw R in any object. For most sources (6/8 short time 
scale and 6/7 long time scale observations), there was sig- 
nificant (>99.7% confidence) evidence of the need to thaw 
r. There was also significant (>99.7% confidence) evi- 
dence of the need to thaw ¥ Ka in the short time scale ob- 
servations of MCG 6-30-15 (I), NGC 3516 and NGC 7469 
as well as in the long time scale observations of MCG 6-30- 
15 (II), NGC 3516 and NGC 4151. These results initially 
suggest statistically significant variations in the line in at 
least some of these objects. Table 5 lists the mean values 
of the 2-10 keV fiux F2-1Q, F and Vxa derived from the 
time-resolved spectral fitting. F2-W was measured from 
the model fits. 

3.3. Error Estimation 
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For the majority of the individual time-resolved fits, the 
resulting values of the fit statistic xl were rather low {xl 
< 1.0 for 243 of the 316 total time bins). As discussed 
by Nandra et al. (2000) and VEOl, this was most likely 
due to the standard software's overestimation of the er- 
rors in any given energy bin in the net spectrum. These 
errors are produced by combining the errors in the source 
with those of the modelled background. The software uses 
Poisson statistics, estimated from a large quantity of back- 
ground data, to estimate the background error. For short 
exposures, this results in an overestimate of background 
errors, which should typically be negligible compared to 
the error on the total (source plus background) rate. To 
better ascertain the variability of the measured spectral 
parameters, the standard errors produced by the software 
were discarded and new la errors were estimated based on 
the variability properties of the derived light curves. The 
point-to-point variance was used; this is defined as 



continuum and line light curves; 
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where Xi are each of the N total points in the light curve. 
VEOl used simulations to validate the accuracy of these 
errors. These are conservative error estimates; only trends 
in the light curves sustained over several time bins will 
exceed the expected error using this formulation. These 
errors are listed in Table 5 for F and Fxa and used here- 
after. Errors on values of F2-10 within a time bin were 
determined from the 2-10 keV continuum light curve, us- 
ing the mean flux error on the ~1 ks exposures in that 
time bin. 

3.4. Time-Resolved Broadband Spectral Fits 

Figures 1 and 2 show light curves of F2-1Q, F, and 
Fxa for the single-Gaussian model on short and long 
time scales, respectively. The results of the Laor model 
time-resolved fits showed identical variability patterns 
and nearly identical time-averaged values to the single- 
Gaussian fits, and they are not plotted. 

This spectral fitting was performed keeping the value of 
R fixed, i.e. with the normalization of the PEXRAV reflec- 
tion component varying in sync with the normalization of 
the power-law. This was done as there was no formal sta- 
tistical requirement to allow the value of R to vary in the 
fits; freezing R also minimized the infiuence of this free 
parameter on measurement of T and Fko- However, as 
the reflection continuum and Fe line are likely to be 
produced in the same location (e.g., the accretion disk) it 
would perhaps be expected that the reflection component 
tracks variations in the line rather than the power-law con- 
tinuum. Time-resolved spectral fitting was repeated with 
the normalization of the reflector tied to the Fna instead. 
The resulting light curves of F^ „ were essentially the same 
as those with a fixed value of R. 

3.5. Fractional Variability Amplitudes 

Fractional variability amplitudes [F-uar] e.g., Vaughan et 
al. 2003, Edelson et al. 2002) were measured to quantify 
the instrinsic variability amplitude relative to the mean 
count rate and in excess of the measurement noise in the 
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where 5^ is the total variance of the light curve, {crl^.^) is 
the mean square error and {X) is the mean count rate of 
A'^ total points. The error on Fyar is 
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This error formu- 



as disciissed in Vaughan et al. 
lation estimates ctf^^^ based on random errors in the data 
itself, and not due to random variations associated with 
red-noise processes. Furthermore, these estimates arc only 
an approximation to the true values of Ui?^^^ since the 
data errors used here instead have been estimated from 
the point-to-point variance. Table 6 lists the continuum 
and line variability amplitudes. 

3.6. Cross Correlation Analysis 

Temporal cross-correlation functions were calculated for 
F2-io~FAraj F2_io-F, and T-Fna using both the discrete 
correlation Function DCF; Edelson & Krolik 1988) and the 
interpolated correlation function (IGF, White fc Peterson 
1994). IGF lag uncertainties were calculated using the 
bootstrap method of Peterson et al. (1998); however, this 
method often underestimates the true errors introduced 
by uneven sampling, e.g., Edelson et al. (2001). Figures 
3 and 4 and Tables 7-9 show the cross-correlation results, 
where a positive shift denotes the first quantity leading 
the second. The DGF and IGF time bins used were 1.0 
and 0.5 times the time bin size in the time-resolved fits, 
respectively. 

One must be skeptical of cross-correlation analysis in 
these circumstances, however, due to the red-noise nature 
of the continuum and line fiuxes (as speculated below). 
Gross-correlation between two unrelated red-noise noise 
light curves can randomly yield spurious correlations with 
higher than expected values of the correlation strength 
Tmax (e.g., Welsh 1999). For instance, the rmax = 0.7 
optical-to-X-ray lag in NGG 3516 claimed by Maoz et al. 
(2000) was not confirmed with additional data by Maoz et 
al. (2002). Gombincd with the fact that most of these light 
curves only have 10 20 points, spurious correlations and 
anticorrelations, including those based on very few points 
(i.e., trends much shorter than the total duration), can 
arise. Discarding peak lags greater than ~l/3 of the total 
duration can eliminate some of these spurious lags; this 
has been done for the values reported in Tables 7-9 and 
the functions plotted in Figures 3 and 4. However, as has 
been done in the present work, it is more practical to not 
trust any lag where rmax is not very close to 1.0, applying 
relatively more skepticism for lower values of rmax- 

Gross-correlation analysis of F2_io and Fxa in the short 
time scale observations yields peak correlation coefficients 
in the range 0.17 0.72. These values do not readily con- 
vince us that the analysis has yielded 'believable' corre- 
lations. For instance, the highest correlation c;ocfficient, 
fmax = 0.72, arises in NGG 3783, which only has 12 
points and for which substantial variability in the line is 
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not detected (see Table 4). NGC 5548, for which cross- 
correlation analysis yields the second-highest value of Tmax 
in the short time scale sample, is a similar situation: the 
reported lag is far from zero, based on only a handful of 
points, and one cannot conclude that this is a reliable cor- 
relation robustly detected in excess of the long time scale 
red-noise variations present. It follows that the lag delays 
and lag errors listed in Table 7 must also be regarded with 
a large dose of skepticism. 

For long time scale Ka correlations, the values 

of Tmax lie in the range 0.45—0.87. Despite the presence 
of higher values of Tmax-, however, caution is again war- 
ranted. For instance, as the continuum and line light 
curves are both dominated by large decreases in NGC 3516 
and NGC 5548, the resulting cross-correlation functions 
have moderately-high values of the correlation coefficient 
at nearly all lags tested, and the peak correlation coeffi- 
cients do not provide evidence of a strong correlation "in 
excess" of this artefact. If there were an increase in the 
values of Vmax towards a particular time scale, this would 
indicate more of a physical relevance associated with that 
time scale. However, only three objects show an increase 
in their respective values of average DCF and ICF Tmax 
towards longer time scales. 

T-F/fa cross-correlation functions show similar correla- 
tion coefficient values as the Y^-w-PKa cross-correlation 
functions. However, many cross-correlation functions be- 
tween F2-10 and r on both time scales yield higher values 
of Tmax- The values of rmax on short time scales span 0.54- 
0.92; those on long time scales span 0.31-0.94 (0.62-0.94 
if NGC 4151 is excluded). All but one of the lag delays 
are consistent with zero. 

3.7. Fyar Spectra 

Fvar spectra, measuring F^ar of the net (background- 
subtracted) count rate in each channel, were calculated as 
a relatively less model-dependent assessment of the broad- 
band spectral variability. Such analysis allows one to study 
the behavior of the varying components. Source and back- 
ground spectra were derived over the time intervals listed 
in Table 10 with the goal of constructing light curves of the 
flux in each channel for at least ~20 points. In addition 
to the observations used in the above analysis, F^ar spec- 
tra were measured for an additional four objects on short 
time scales and five objects on long time scales; the total 
of twelve observations on each time scale included nine ob- 
jects common to both lists. Because the PCU gain settings 
changed three times since launch and the channel energy 
boundaries differ between gain epochs, the longest dura- 
tion that did not cross a gain epoch boundary was used. 
The error on each light curve point was assigned using the 
point-to-point variance (see §3.3). The highest-energy bins 
were re-grouped as necessary to insure adequate signal-to- 
noise in all energy bins. For the soft-spectrum (F ~2.7; 
Vaughan ct al. 1999) source Akn 564, data above 9 kcV 
were discarded to minimize the influence of systematic 
variability in the RXTE background models. Similarly, in 
Fairall 9 (F - 2.2 over 2-10 keV, TARTARUS database), 
data above 14.3 keV were discarded for the short time 
scale data. F^ar and (Jf^^t were calculated in each chan- 
nel using the equations in §3.5. The final F^ar spectra 
are presented in Figures 5 and 6 for short and long time 



scales respectively. For the seven sources with published 
power spectral densities (Markowitz et al. 2003, Uttley, 
McHardy & Papadakis 2002), the relative continuum vari- 
ability amplitudes obtained on short and long time scales 
are in reasonable agreement with values predicted from 
their published power spectral density shapes. 

4. RESULTS 

A summary of the observed continuum broadband flux 
and spectral variability properties of the sample is given in 
§4.1. The observed line variability properties are summa- 
rized in §4.2. Finally, a comparison between the continuum 

and line variability properties is given in §4.3. Notes on 
individual sources are deferred to Appendix A. 

4.1. Continuum variability 

All objects display significant variability in continuum 
flux: in x^-tests against the constant hypothesis, the con- 
stant hypothesis is rejected at >99.99% confidence for all 
continuum light curves. Visual inspection of Figures 1 
and 2 confirms this. There is evidence for the continuum 
to vary more strongly towards longer time scales: of the 
six objects with both short- and long time scale monitor- 
ing, the value of continuum F^ar increases towards longer 
time scales for five objects. The F^ar spectra also ex- 
hibit stronger continuum variability on longer time scales 
in most sources. Both these analyses support an increase 
in variability amplitude towards longer time scales that 
is consistent with the long-established red-noise nature of 
Seyfert 1 continuum X-ray emission. 

The majority of F^ar spectra show the softer energies to 
be the most variable, consistent with softening of the dom- 
inating broadband X-ray continuum as sources brighten 
(e.g., Markowitz & Edelson 2001, Nandra et al. 1997). The 
only exception is Akn 564, which exhibits a much flatter 
dependence on F^ar with energy. The presence of F2_io- 
F correlations with no measureable lags in many objects 
on both time scales supports this idea. Such behavior has 
been observed numerous times before (e.g., VEOl, Done 
et al. 2000, Lee et al. 2000, Chiang et al. 1999, Nandra, 
Pounds & Stewart 1990). 

Above ^10 keV, unmodelled systematic variations in 
the background can contribute greatly to the total ob- 
served variability. In most sources, therefore, the influence 
of the variability of the Compton reflection hump on the 
broadband continuum above ~10 keV cannot be reliably 
constrained from the Fyar spectra. 

4.2. Fe Ka line variability 

Given that the 5-7 keV continuum flux greatly outshines 
the line flux and given the nature of spectral fitting, it is 
possible that there could exist some systematic effect that 
corrupts the determination of the continuum level at the 
~few per cent level, in turn corrupting determination of 
the line fiux by ^10 percent or higher. However, such 
an effect is not readily identifiable, and the line variations 
measured are hence assumed intrinsic to the sources. 

Timc-rcsolved spectral fitting has revealed a wide range 
of variability in the line on both time scales. On short time 
scales, MCG-6-30-15 (I) shows evidence for line variabil- 
ity: a x^-test against a constant hypothesis results in re- 
jection of that hypothesis at >99% confidence. Similarly, 
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the constant hypothesis is rejected in MCG-6-30-15 (II) 
and NGC 3516, but only at 98% and 95% confidence, re- 
spectively. 

On long time scales, the number of objects showing line 
variability increases: in three of the seven observations 
(NGC 3516, NGC 4151 and NGC 5548), the constant line 
flux hypothesis is rejected at >99.99% confidence. In two 
additional observations, MCG-6-30-15 (I) and 3C 390.3, 
the constant line flux hypothesis is rejected, but only at 
97.9% and 97.7% confidence, respectively. The most note- 
able trends in line flux are the gradual decreasing trends 
throughout the long time scale observations of NGC 3516, 
NGC 4151 and NGC 5548. Lack of Hne variability cannot 
be ruled out in MCG-6-30-15 (II) and NGC 3783. Over- 
all, one can conclude that Fe Kq, line variability is more 
commonly detected towards longer time scales. 

This tendency is conflrmed using measured values of 
F-uar ■ the value of line Fyar increases towards longer time 
scales for four objects, the values of line Fyar on both time 
scales are consistent for the remaining two objects, and in 
no case does line Fyar decrease towards longer time scales. 
It is therefore plausible to conclude that, at least on the 
time scales probed herein, the iron line variability is, like 
the continuum variability, a red-noise process. 

The average maximum peak-to-trough line flux varia- 
tion detected is a factor of ~2 for all observations; typi- 
cally, this is a factor of ~4.5 greater than the typical er- 
ror bar. The largest peak-to-trough variations 
in the short-term observation MCG-6-30-15 (II) and the 
long-term observations of NGC 4151 and NGC 5548. The 
smallest peak-to-trough variations are ^40-50%, as with 
the short-term observations of NGC 3516 and NGC 3783 
and the long-term obsrvations MCG-6-30-15 (II) and 
NGC 3783; typically this is a factor of ~2 greater than 
the typical error bar. With the time bins chosen in this 
experiment to maximize the variability-to-noise in the iron 
line light curves, this experiment may not be sensitive 
to smaller amplitudes of variability (peak-to-trough varia- 
tions of ~30% or less). 

4.3. Relating continuum and line variability 

Visual examination of the continuum and line light 
curves does not reveal any obvious relation between them 
in most objects. Large-amplitude trends in continuum of- 
ten go unmatched by simultaneous changes in the line; 
prominent examples of this are listed for individual sources 
below. In only a few observations do the continuum and 
line appear to track each other: in the long time scale 
observation of NGC 5548, both the continuum and line 
decrease by a factor of ^3 over the 4-year duration. There 
are also somewhat similar decreasing trends dominating 
both the continuum and line in the long-time scale obser- 
vations of NGC 3516 and NGC 4151. 

Measured variability amplitudes, Fyar spectra, and the 
increased detection of line variability in the light curves to- 
wards longer time scales all indicate that both the contin- 
uum and line appear to be more strongly variable towards 
longer time scales in most objects. However, there is evi- 
dence that the line does not vary as strongly as the contin- 
uum. In 11/15 observations, the value of continuum Fyar 
is greater than the line Fyar- For the eight of these eleven 
observations with a defined value of line Fyar (positive ex- 



cess variance) , the continuum and line excess variances are 
compared via an F-test. In six of these eight observations, 
the null hypothesis of the continuum and line excess vari- 
ances being consistent is rejected at >98.2% confidence. 
In four of the fifteen observations, the values of continuum 
and line Fyar are consistent with each other. For three 
of these observations, an F-test cannot reject the null hy- 
pothesis of the continuum and line excess variances being 
consistent at greater than 79% confidence; it is rejected at 
>99.99% confidence in the fourth observation (the long- 
term NGC 3783 observation). There are no observations 
with a larger value of Fyar for the line compared to the 
continuum. 

As seen in Figures 5 & 6, in many of the Fyar spectra, 
there is a reduction in variability at ~6 keV. The average 
drop in Fyar at ~6 keV relative to the surrounding contin- 
uum is ~10%, or ~5 times larger than the typical errors. 
Best fit lines to the Fyar spectra including and excluding 
the ~5-7 keV region indicate that the dip is significant 
at greater than 90% confidence in an F-test for 11/24 ob- 
servations. These dips are consistent with dilution of the 
~5-7 keV continuum variability by relatively reduced line 
variability. Such a feature has been noted previously in 
the Fyar spectra of NGC 4151 (Schurch & Warwick 2002) 
and MCG-6-30-15 (Fabian et al. 2002). The presence of 
this feature supports the argument for decoupling of the 
line and continuum variability properties. Additionally, 
cross-correlation analysis does not reveal any significant 
temporal correlations between continuum and line varia- 
tions. 

Overall, there hence is little support for a simple or di- 
rect temporal relation between the observed continuum 
and Fe Kq, line variability properties. This result broadly 
agrees with the study of Weaver, Gelbord & Yaqoob 
(2001), who did not find any evidence that Fe K„ line 
variability on time scales of days to years was correlated 
with continuum variations; their ASCA sample inchided 
all of the objects in this paper's time-resolved spectroscopy 
samples except 3C 390.3, though with much less continu- 
ous sampling. 

5. DISCUSSION 

This paper examines the variability properties of both 
the continuum and Fe Kq, line on time scales from days 
to years. The continuum variability properties observed 
in the current sample have been observed before in many 
Seyfert 1 galaxies. The data show an increase in broad- 
band variability amplitude towards longer time scales, 
demonstrating the red- noise nature of the continuum vari- 
ability. Additionally, there is an increase in variability am- 
plitude towards the softer X-rays, consistent with sources' 
softening as they brighten. In contrast, the Fe Kq, line 
shows a divergence of behaviors. Most observations show 
no clear trend whatsoever. Only a few observations show 
temporal trends in the line fiux that mimic the contin- 
uum to at least some degree, and only on time scales of 
years. As discussed in §4.3, however, the line varies less 
strongly than the continuum in many sources, as seen from 
a comparison of the spectral fit line and continuum flux 
light curve variability amplitudes as well as from a reduc- 
tion in broadband Fyar near 6 keV seen in many sources' 
Fvar spectra. The relationship between continuum and 
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line variability is examined further in §5.1. The broad- 
band spectral variability results are discussed in §5.2. 

5.1. Fe Ka line variability 

This paper has demonstrated the presence of a range of 
variability behavior in the Fe Kq, line. The results permit 
tests of simple models in which the bulk of the continuum 
originates in a central corona and the line variations are 
driven solely by continuum flux fluctuations modiflcd by 
light-travel time effects. For instance, such models imply 
that there should exist a connection between the origin of 
the line emission (which can be inferred from the resolved 
profile shape) and the presence of rapid line variability. If 
the bulk of the line is emitted in the inner accretion disk 
(as implied by the broad line profiles in MCG-6-30-15 and 
NGC 3516), then the continuum and line variations should 
be correlated on time scales comparable to the light cross- 
ing time of the inner accretion disk. This is because the 
reprocessing times arc likely very short compared to the 
light travel times, which dominate the correlation. How- 
ever, the time bins used in this paper are typically at least 
two orders of magnitude longer than the expected light 
crossing time scale of the putative inner accretion disk. 
This means that any variations on this time scale would 
be grossly undersampled. 

If the bulk of the line is emitted far from the central 
corona (e.g., as implied by the narrow line profiles for 
NGC 4151 and NGC 5548), then one would expect the 
line variations to track the continuum variations only on 
sufficiently long time scales (e.g., longer than days-weeks 
if the line is emitted in a region coincident with the broad 
line region; Kaspi ct al. 2000). For instance, the presence 
of similar long-term trends in the continuum and line in 
NGC 5548 supports a model in which the bulk of the line 
emission originates far from the central corona in this ob- 
ject. More specifically, it is consistent with the prediction 
of Yaqoob ct al. (2001) that the flux of the narrow core in 
NGC 5548, estimated to originate in a region consistent 
with the broad line region in this object, should track the 
continuum flux on time scales of ^5 - ^100 days. How- 
ever, such trends are not ubiquitous. For Fe Kq, lines orig- 
inating far from the central corona, one would expect to 
see an increase in number or strength of continuum-line 
correlations as one goes towards towards long time scales. 
However, there is no evidence from the current sample to 
support this. 

The lack of correlation between the continuum and the 
line suggests that modification to this straightforward pic- 
ture is required. One possible explanation for the reduced 
variability in the line relative to the continuum is that a 
portion of line emission originating in the inner disk is ab- 
sent or suppressed. The inner portions of the accretion 
disk may be truncated, highly ionized, or radiatively inef- 
ficient (e.g., as in an advection-dominated accretion flow; 
Narayan & Yi 1994). Geometry may play a role: a rel- 
atively more concave disk, for instance, can reprocess a 
larger fraction of continuum photons in its outer regions 
(Blackman 1999). For cases in which the bulk of the ob- 
served line flux originates in the outer disk, the direct re- 
sponse of the line peak to rapid continuum variations is 
reduced. Another possibility is that if the X-ray contin- 
uum emission is produced in a corotating flare some height 



above the reflecting disk, then the orbital motions of the 
flare and disk can also lead to complex temporal behav- 
ior in the observed line flux (Ruszkowski 2000). Fabian 
& Vaughan (2003) suggested that light-bending near the 
black hole can lead to substantial de-coupling of the ob- 
served continuum and line fluxes. This happens because 
much of the emission from the X-ray source is focussed on 
to the disc whereas only a small fraction is able to reach the 
observer. In such a scenario small changes in the geometry 
of the X-ray source can lead to relatively large changes in 
the amount of continuum radiation reaching the observer 
(hence the strong continuum variations) but relatively lit- 
tle change in the total amount of radiation illuminating 
the disk. 

Another intriguing possibility relies on the presence of 
an ionized skin on the disk surface. Thermal instabilities 
can lead to formation of such an ionized skin (Ross, Fabian 
& Yoimg 1999; Nayakshin, Kazanas & Kallman 2000) with 
most of the observed neutral line emission originating in 
the colder layers beneath. The ionized skin can scat- 
ter incoming continuum photons and outgoing line pho- 
tons, affecting the observed neutral line flux, as discussed 
by Nayakshin (2000) and Nayakshin & Kallman (2001). 
However, changes in the ionization state of the skin are 
linked to the dynamic time scale for readjustment of hy- 
drostatic equilibrium (Nayakshin & Kazanas 2002 for a 
lamppost illuminating model; see also Collin et al. 2003 
for a multiple-flare illuminating model). This time scale 
is usually an order of magnitude or more larger than the 
light- crossing time; this can be >minutes-hours if most of 
the line emission originates in the inner regions of a typ- 
ical thin accretion disk (e.g., around a 10®~* Mq black 
hole). The observed line flux hence may not respond to 
continuum variations shorter than this time scale, result- 
ing in decorrelation of the observed continuum and line 
flux variations. Furthermore, the ionization statc^ of tlit^ 
skin depends on the high-energy cutoff of the power-law 
continuum (Nayakshin & Kallman 2001); this could be 
variable, but it is difficult to constrain with RXTE. Inter- 
estingly enough, such a variation in this parameter could 
lead to a change in Fe K^ line flux even without a change 
in 2-10 keV continuum flux. However, recent studies have 
not been able to directly link the line flux behavior with 
the ionization state of the disk, e.g., Ballantyne, Vaughan 
& Fabian (2003). 

All of the above scenarios are still consistent with the 
prediction that both broad and narrow Fe Kq, lines should 
eventually track continuum variations on time scales of 
months to years. The finding that the line decreased along 
with the continuum in the long time scale observation of 
NGC 5548 is consistent with this idea; however, such cor- 
related behavior is not seen throughout the entire sample. 
For line emission originating in the outer disk or molecu- 
lar torus, ionization in the reflecting material is negligible 
due to the low X-ray continuum flux. However, other pro- 
cesses may be relevant for decorrelating long-term contin- 
uum and line variations. For instance, it may be possible 
for an ionized wind (e.g., George et al. 1998) near the site 
of reflection to mimic an ionized 'skin' in these regions. 
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Most of these objects show a strong correlation with 
no measureable lag between 2-10 keV continuum flux and 
photon index. Additionally, Fyar spectra show that in 
nearly all the objects, the soft-spectrum (narrow line) 
source Akn 564 being an exception, variability amplitudes 
tend to increase towards relatively softer energies. These 
observations show that the broadband emission steepens 
as it brightens, in line with numerous similar observations 
(e.g., Magdziarz et al. 1998, Done et al. 2000). 

The tinic-resolvcd spectral fitting of MCG 6-30-15 per- 
formed from ASCA data by Shih, Iwasawa & Fabian 

(2002) and from XMM-Newton data by Fabian & Vaughan 

(2003) point towards a scenario in which flux-correlated 
changes in the overall 2-10 keV continuum slope are ex- 
plained by the superposition of a soft, variable power-law 
likely associated with coronal emission and a harder, much 
less variable spectral component likely associated with the 
Compton reflection hump. Taylor, Uttlcy & McHardy 
(2003) use flux-flux plots to show that such a description of 
spectral variability works in MCG-6-30-15 and NGC 3516. 
Specifically, in this model, the intrinsic slope of the under- 
lying power-law does not vary significantly; rather, it is 
the relative contributions of the two spectral components 
that change the overall spectral shape as observed. This 
model does not require physical properties of the power- 
law emitting corona such as electron temperature, optical 
depth, or flare size to vary signiflcantly. Furthermore, this 
two-component model may have relevance to explaining 
the observed Fe variability. As the Fe line and 
Compton hump are expected to originate together, a rel- 
atively constant Compton hump could be associated with 
variability in the Fe Kq, line that is reduced relative to the 
continuum variability, as observed (e.g., Shih et al. 2002). 

Alternate models to explain the observed spectral vari- 
ability rely on flux-correlated changes in the power-law 
slope of the coronal component. In the context of popular 
thermal Comptonization corona models (Haardt, Maraschi 
& Ghisellini 1997, Poutanen & Fabian 1999), a flux- 
photon index correlation can arise if the X-ray variability 
is due to changes in the input seed photon population. In 
contrast to the two-component model above, some phys- 
ical property of the corona may vary with the 2-10 keV 
X-ray flux in these types of models. 

The Fyar spectra of the soft-spectrum (narrow-line) 
Seyfert 1 galaxy Akn 564 display a much flatter depen- 
dence on energy in the 3.5 9 keV bandpass (and may 
even be slightly increasing with energy on the long time 
scales) relative to the other objects, most of which are 
hard-spectrum (broad- line) Seyfert 1 galaxies. This has 
been seen previously in Akn 564 by Edelson et al. (2002) 
as well as in other soft-spectrum sources (Vaughan et 
al. 2002). This could imply that, in the context of the 
spectral-pivoting model above, the broadband spectral 
shape changes very little in soft-spectrum sources, per- 
haps due to the lack of variation of some physical pa- 
rameter that does vary in broad-line Seyfert 1 galaxies. 
Alternatively, if the two-component model above is rele- 
vant, it is possible that soft-spectrum sources have much 
weaker hard/constant components, leading to relatively 
little change in the overall spectral shape with flux. Over- 
all, however, it is difficult to present a simple, unified phe- 
monological model to simultaneously explain the depen- 



dence on energy of the variability in broad-line Seyfert 1 
galaxies and lack thereof in narrow-line Seyfert 1 galaxies. 

6. CONCLUSIONS 

Time-resolved broadband spectral fitting has been sys- 
tematically performed on near-continuous RXTE monitor- 
ing data of seven Seyfert 1 galaxies in order to probe the 
behavior of emission components on time scales of both 
~days to ~weeks and ~months to ~years. Fe K^j line 
variability is present at a variety of levels. It is readily 
detectable in some observations; for instance, the average 
peak-to-trough variations are a factor of ~2 and as high 
as ~3 in three cases (the short-term observation MCG- 
6-30-15 (II) and the long-term observations of NGC 4151 
and NGC 5548). In other observations, the line does not 
display appreciably variability: pcak-to-trougli variations 
are only ~40-50%. The measured values of Fyar indicate 
that a higher fraction of the objects sampled show de- 
tectable line variability towards longer time scales. Of par- 
ticular interest are systematic long-term decreases in line 
fiux over durations of 3-4 years in NGC 3516, NGC 4151, 
and most noteably NGC 5548. Furthermore, the fractional 
variability amplitude of the line increases towards longer 
time scales in many objects, "tracking" the red-noise con- 
tinuum. However, fractional variability amplitudes show 
that the line does not vary as strongly as the continuum. 
Fyar spectra also generally show a dilution in ~5-7 keV 
continuum variability due to the relatively less variable 
Fe Kq, line. In NGC 5548, the continuum and line simul- 
taneously decrease by similar factors over a 4-year period; 
the long-term observations of NGC 3516 and NGC 4151 
also show somewhat similar trends in continuum and line 
flux. For the sample as a whole, however, there appears 
to be no simple temporal relation between the continuum 
and the line, contradicting simple models in which the 
line flux is driven solely by continuum variations. Local 
effects such as an ionized skin may be relevant in decorrc- 
lating the continuum and line variations, particularly for 
lines originating in the inner accretion disk. The strongest 
statement that can be made on the basis of these results is 
that for most of the Seyfert 1 galaxies studied, the Fe K„ 
emission is less variable than the continuum, and in those 
cases where there does appear to be evidence for changes 
in the line flux, these changes are not well correlated with 
the continuum. 

There is evidence for a correlation with no measure- 
able lag between continuum flux and measured photon in- 
dex on both time scales in many sources. Relatively less 
model-dependent Fyar spectra generally show the broad- 
band continuua to be more variable at relatively softer 
energies, supporting Seyfert 1 galaxies' softening as they 
brighten, though the soft-spectrum source Akn 564, with 
its flat Fyar spectra, is an exception. 
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APPENDIX 

NOTES ON INDIVIDUAL SOURCES 

• 3C 390.3: The short time scale timc-rcsolvcd and time-average fit results arc in reasonable agreement with Gliozzi ct 
al. (2003): the light curve for P has absolute values slightly systematically lower by <0.1 than Gliozzi et al. (2003) but 
the variability patterns are identical. On short time scales, a decrease in F2_io by a factor of ~2 over a span of 30 days 
does not appear to have any simultaneous effect on F/fai which remains nearly constant during the observation. On long 
time scales, there is a doubling in F2-10 over 300 days, and an increase of ~50% in Fkc over approximately the same 
period. 

• MCG-6-30-15: Due to improved background models, the current analysis yields a lower absolute value for R for 
the integrated spectrum of the short time scale observation MCG-6-30-15 (I) than in VEOl (values of R obtained from 
spectral fitting arc highly sensitive to the shape of the modelled background spectrum above ^10 keV), but the time- 
resolved spectral analysis and the cross-correlation analysis on the short time scale observation MCG-6-30-15 (I) is in 
excellent agreement with VEOl. Fku shows variability on short time scales, as well as in the long time scale observation 
MCG-6-30-15 (I). The long time scale observation MCG-6-30-15 (II) does not yield a significant detection of variability in 
the line; one possibility is that on these time scales probed (~100 d to ~2000 d) the variability is "averaged out" to some 
degree. There arc examples in these observations of large increases or decreases in F2-10 (in some cases by a factor of ~2) 
that do not have any obvious immediate effect on Fxa- The 6.4 keV dips in the Fyar spectra are prominent, especially for 
the short time scale observation MCG-6-30-15 (I). The four observations, when considered in order of increasing duration, 
do not show any obvious pattern to the Fe line variability. 

• NGC 3516: On short time scales, F2_io and Fxa exhibit uncorrelated variability. Persistent trends in F2_io such as 
the rapid decrease from Modified Julian Day (MJD) 50550-50554 are not simultaneously matched by similar behavior in 
F/fa- Rapid changes in F/fa occur without sudden changes in F2_io, e.g., the decrease in F^a by ~20% at MJD 50551 
and also the increase in F^a by ~30% at MJD 50558 as illustrated in Figure 7. The 6.4 keV dip in the short time scale 
Fyar spectra is especially noticeable. On long time scales, both F2_io and Fxa show an overall decrease by a factor of 
~2 over the 1000-day duration. 

• NGC 3783: Formally, F^a does not demonstrate appreciable variability on either time scale. Relatively rapid changes 
in F2_io by >10% towards the middle of the short time scale observation and ~20% - ~40% towards the beginning and 
end of the long time scale observation do not appear to have any major simultaneous impact on F^a ■ 

• NGC 4151: On short time scales, Fxa is not formally variable. One note of interest, for instance, is the doubling of 
F2-10 during the last 4 days of the short time scale observation while Fxa remains virtually constant during that time. 
The 6.4 keV dip in the short time scale Fyar spectra is especially noticeable. On the long time scale, Fku is much more 
variable. There are similar decreasing trends by a factor of ~3 in both the continuum and line over the whole 800-day 
duration. Of note is the sudden dip and rise in the line around MJD 51620 in the long time scale observation while F2_io 
remains nearly constant; Figure 7 demonstrates the dip. 

• NGC 5548: On long time scales, a gradual decreasing trend by a factor of 3-4 dominates the behavior of both F2_io 
and Fkc- The Fe line variability amplitude is substantially larger on the long time scale than the short due to this 
dominating long time scale trend. 

• NGC 7469: All the time-average and time-resolved spectral fit results and cross-correlation results are in reasonable 
agreement with Nandra et al. (2000). The current analysis yields absolute values of P that are slightly systematically 
lower by <0.1 and the absolute values of Fxa obtained in the current analysis are systematically ^^30% higher compared 
to Nandra et al. (2000), but the variability patterns are consistent. 
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Table 1 

Source and sampling parameters 



Source 


2-10 keV Lum. 




MJD Date 


Mean 


Spectral Fits 


Spectral Fits 


Name 


(log(erg s"^)) 


z 


Range 


c s~-^ 


Time Bin 


No. Pts 


3C 390.3 


44.26 


0.056 


50220.6 50278.0 


2.3 


6.0 d 


10 


MCG 6-30-15 I 


42.80 


0.008 


50664.1 50672.5 


4.1 


0.13 d 


63 


MCG 6-30-15 II 


42.83 


0.008 


51622.7 51676.0 


4.4 


1.6 d 


34 


NGC 3516 


42.90 


0.009 


50523.0-50657.4 


4.2 


6.4 d 


21 


NGC 3783 


43.12 


0.010 


51960.1-51980.0 


5.7 


1.7 d 


12 


NGC 4151 


42.24 


0.003 


51870.6-51904.9 


7.4 


1.4 d 


24 


NGC 5548 


43.52 


0.017 


52091.6 52126.7 


4.5 


3.2 d 


11 


NGC 7469 


43.2.J 


O.OIG 


50244.0 50275.7 


2.8 


1.3 d 


25 


3C 390.3 


44.40 


0.056 


51186.0 51964.8 


3.2 


51.9 d 


15 


MCG 6-30-15 I 


42.81 


0.008 


51926.7-52252.9 


4.2 


23.4 d 


13 


MCG-6-30-15 II 


42.85 


0.008 


50159.8-52252.9 


4.6 


128.7 d 


16 


NGC 3516 


42.86 


0.009 


50523.0-51593.4 


3.8 


53.5 d 


20 


NGC 3783 


43.22 


0.010 


51180.5-52375.2 


7.1 


59.8 d 


20 


NGC 4151 


42.55 


0.003 


51179.6-51964.6 


14.9 


43.6 d 


18 


NGC 5548 


43.55 


0.017 


51183.9-52643.1 


4.8 


104.2 d 


14 



Note. — The short time scale data arc listed in the top half of the tabic; the bottom half hsts the long 
time scale data. Column (2), log of the 2-10 keV luminosity, was calculated using the online W3PIMMS tool 
and using the mean 2 10 kcV count rate per PCU, X-ray photon indices and neutral absorption based on the 
time-averaged spectral fitting below, and assuming Hq = 75 km s~^ Mpc~^ and qo = 0.5. The count rates in 
Column (5) are 2-10 keV count rates per PCU. 
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Table 2 

Time-average fit parameters, single-Gaussian model 



Source 






Eo 


cr 








Name 


Nh 


r 


(keV) 


(keV) 


FKcr 


R 


xVd.oi. 



3C 390.3 4.28 

MCG-6-30-15 1 4.09 

MCG-6-30-15 II 4.09 

NGC 3516 2.94 

NGC 3783 8.26 

NGC 4151 616 

NGC 5548 1.75 

NGC 7469 4.86 



1.641±0.006 
2.018±0.010 

1.947±0.016 
1.573±0.007 
1 666+°°22 

i.DOD_Q gig 

1.524±0.010 
1.733±0.017 
1.845±0.009 



6 00+° °'* 
5.87±0.02 



6.17 



+0.02 
-0.03 



6.02±0.01 
fi 27+° °'* 

"•^ ' -0.03 

6.40±0.01 
6.27±0.06 
6.31±0.02 



0.44±0.07 
0.68±0.02 

0.46±0.04 
5Q+o °i 

0.41±0.06 

0.42±0.10 
0.28±0.04 



0.64±0.06 

1.79±0.06 

1.35±0.07 
9 (qo+o.oe 

^•"^-0.04 

1.75±0.12 

3.45±0.08 
n 70+0.07 

'-'•'''-0.08 

0.64±0.02 



< 0.03 
0.85±0.07 
1.15±0.12 
20+° °'* 
0.62±0.12 
0.55±0.05 
0.44±0.10 
0.42±0.05 



1.0+0_g 

2 020+° °3° 

Z.UZU_o.o29 

2.000±0.006 

l.D+Z_o.o04 
1 RQQ+0.008 

J^.oyo_o 007 
1 692+°°" 

i.ut;z_0 005 

1.739±0.010 



-+0.O8 

-•-0.11 
hO.lO 
-0.09 
hO.Ol 
="-0.02 
h0.02 
-0.01 



-O.Ofj 
-■-0.02 
-+0.28 
-0.25 
-1 17+0.05 
J---'- '-0.04 

0.27±0.03 

7+0.05 
'-0.04 

fS5+°-°'* 



61.0/34 
292.0/33 

23.0/27 
487.6/33 

22.3/28 
101.0/26 

25.3/28 
146.0/33 



3C 390.3 4.28 

MCG-6-30-15 1 4.09 

MCG-6-30-15 II 4.09 

NGC 3516 2.94 

NGC 3783 8.26 

NGC 4151 812 

NGC 5548 1.75 



6.10±0.04 
t- 02+0 07 

6.00±0.01 

6.09±0.01 
6.23±0.02 
6.40t°o;°? 
6.28t°;°l 



0.25^ 

0.75^ 

0.59j 

0.47^ 

0.43±0.03 

i2+° °2 
'-'■■'-^-0.12 



0.36 



+0.05 
-0.06 



0.47±0.04 
1.76±0.02 
1.75±0.03 
2.11±0.03 
2.04±0.06 

q 79+0.04 
' ^-0.08 
n '7fi-|-0.04 



0.03: 

1.96^ 



0.47] 



0.41 



+0.02 

-0.08 



122.0/97 
30.75/28 
726.9/98 
1009.7/65 
206.6/98 
581.6/94 
105.0/99 



Note. — Results of fits of the model (power-law -|- pexrav -|- Gaussian) x absorption to the time-average 
data for all targets. The short time scale data are listed in the top half of the table; the bottom half lists 
the long time scale data. The PEXRAV component was omitted for the short time scale spectrum of 3C 390.3. 
Column (2) lists the neutral absorption column in units of 10^° cm~^, modelled using the Galactic column for 
all sources except NGC 4151. The line energy, E,-,, Column (4), is the energy in the observed, not the rest, 
frame. Column (5) lists cr, the FWHM of the Gaussian. Column (6) lists the Fe Kq, line normalization in 
units of 10^^ ph cm~^ s~*. Column (7) lists i?, the relative reflection of the Compton reflection hump. Errors 
listed above are Icr errors determined by XSPEC. 
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Table 3 

Time- AVERAGE fit parameters, Laor model 



Source 






Eo 












Name 


Nh 


r 


(keV) 


(3 






R 





3C 390.3 4.28 

MCG-6-30-15 1 4.09 

MCG-6-30-15 II 4.09 

NGC 3516 2.94 

NGC 3783 8.26 

NGC 4151 618 

NGC 5548 1.75 

NGC 7469 4.86 



1 RQ7+0.007 
1 QOO+0.019 

i-y^^-o.ois 

1.507±0.004 
1.650±0.022 

1 coA+O.Oll 

1.0Z'i_g QQg 

1 7Oc:+0.017 
1. / ^>J_o.025 



6.19±0.5 

n 47+0.02 

O-'i' -0.03 

6.34±0.3 

n AO + 0.01 

O-4O-0.03 

6.46±0.5 
6.42±0.1 

"•^•J-0.07 



1.832±0.008 6.43±0.2 



l.DOO_o.002 

1.975+"-"2i 



3.2±0.1 
3.6t^o-j 

8.0^}? 
2 5+°-i 

o n+6.9 
''•^-S.O 

2-8^0:8 
2.2±0.1 



13t^ 
2.2i?;« 

13t^ 
7 1+1-1 

' -1-0.1 

c 0+4.8 
O.o_^ 2 

2.8±1.4 



"2fifF™~ 

^DU_200 

8 9+1-" 

1 «+0-4 
^•8-0.5 

9 Q+0.1 

6.3±0.8 

+ 170 

-60 

44^L 



0.68^!];E]^ 
2.76t°;fe 
1.37±0.07 
3.09±0.04 

1-^3-0.37 
o 40+0.10 

n 7q+o:o8 

' ^-003 



^■++-0.03 

1- '"-0.10 
2 64+0-01 

^•"+-0.07 

2 73+0.02 
2.37 



< 0.02 
0.45±0.08 
1.06±0.12 

< 0.03 

0.57lH^ 
0.55±0.05 

40+°-^" 

"•^^-0.09 
j+0.03 
-0.06 
,+0.U6 



0.38"* 



55.6/33 
98.2/32 
19.6/25 
101.5/32 
20.4/27 
101.0/25 
24.5/27 
117.4/32 



3C 390.3 4.28 

MCG-6-30-15 1 4.09 

MCG- 6-30-15 II 4.09 

NGC 3516 2.94 

NGC 3783 8.26 

NGC 4151 690 

NGC 5548 1.75 



1.919 



-0.030 
-0.005 
-0.010 

1.582±0.002 
1 664+"-""^ 

1.00'±_Q QQg 

1 604+°-°^° 

i.DU+_g oil 

1.732±0.011 



6.12±0.4 
6.41±0.9 

6 43+0-01 
"•^"J-0.02 

6.37±0.1 
6 40+0.03 

6.39±0.1 
6.38±0.3 



8.2: 
9.9j 



3.1 

hO.O 

'-6.9 
2.9±0.1 

2.6±0.1 

2 9+0-2 
"'•^-o.i 

8.6l^i 



0.01 

0.22 

-0.07 
-0.03 
-0.05 



230 



-0.09 
0.17 
-0.03 

4.04t°;^ 

0.04 
0.03 



0.06 
1.59 
0.81 



0.08±0.02 

0.38±0.03 

^4+0-02 

"•^+-0.03 



0.73^ 



0.40 



+0.05 
-0.03 



122.1/96 
27.1/27 
457.0/97 
766.8/64 
185.4/97 
450.5/91 
120.5/97 



Note. — Fit parameters for the Laor model fits to the time-averaged data. The short time scale data are listed 
in the top half of the table; the bottom half lists the long time scale data. The pexrav component was omitted for 
the short time scale spectra of 3C 390.3 and NGC 3516. Column (2) lists the neutral absorption column in units 
of 10^° cm~2, modelled using the Galactic column for all sources except NGC 4151. Column (4) is the emissivity 
index (3. Column (5) is the innermost disk radius of line emission, in units of GM/c^ where M is the black 
hole mass. Column (7) lists the Fe line normalization in units of 10~^ ph cm^^ s^i. Column (8) lists R, the 
relative reflection of the Compton reflection hump. Errors listed above are la errors determined by XSPEC. 
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Table 4 

F-TEST RESULTS TO THAW F, Fkc OR R IN THE SINGLE- GAUSSIAN MODEL 



Source 


r 


r 






R 


R 


Name 


F 


Prob 


F 


Prob 


F 


Prob 


3C 390.3 


5.8 


1.2E-10 


1.80 


3.7 E-2 


N/A 


N/A 


MCG-6-30-15 I 


1.78 


1.9E-4 


1.63 


1.4E-3 


0.76 


0.92 


MCG-6-30-15 II 


2.4 


1.6E-5 


1.70 


5.7E-3 


1.31 


0.11 


NGC 3516 


24.5 


<lxlO-''° 


3.20 


2.9E-6 


0.39 


0.99 


NGC 3783 


0.81 


0.64 


0.72 


0.74 


0.53 


0.90 


NGC 4151 


1.94 


4.8E-3 


0.82 


0.71 


1.04 


0.41 


NGC 5548 


1.5 


0.26 


0.99 


0.45 


0.47 


0.92 


NGC 7469 


3.80 


1.7E-9 


2.05 


1.9E-3 


1.58 


3.5E-2 


3C 390.3 


5.84 


3.2E-11 


0.94 


0.52 


0.65 


0.84 


MCG-6-30-15 I 


2.21 


2.9E-2 


1.65 


4.4E-2 


1.80 


2.3E-2 


MCG-6-30-15 II 


2.90 


1.3E-4 


2.29 


2.9E-3 


0.86 


0.62 


NGC 3516 


13.0 


2.3E-36 


4.83 


3.8E-11 


1.08 


0.36 


NGC 3783 


1.17 


0.28 


1.03 


0.43 


0.53 


0.95 


NGC 4151 


12.3 


5.2E-31 


7.74 


2.2E-18 


1.76 


2.7E-2 


NGC 5548 


1.77 


4.11E-2 


1.83 


3.2E-2 


0.85 


0.61 



Note. Results of F-test to determine which parameters to thaw in the 
time-resolved fits. High values of the F-statistic and low values of the proba- 
bility (of observing that value of F from a random set of data) indicate that 
the fits show significant improvement to thaw that parameter. The short time 
scale data are listed in the top half of the table; the bottom half lists the long 
time scale data. 



Table 5 

Mean spectral fit variability parameters and errors 



Source 


Mean F2-10 


Mean 


Mean Fxa 


Name 


(10-" ph cm-2 s-i) 


F 


(10~4 ph cm-2 s-i) 


3C 390.3 


2.58 ± 0.05 


1.633 ± 0.023 


0.639 ± 0.174 


MCG-6-30-15 I 


4.73 ± 0.05 


2.005 ± 0.063 


1.829 ± 0.250 


MCG-6-30-15 II 


4.58 ± 0.10 


1.914 ± 0.103 


1.384 ± 0.261 


NGC 3516 


4.80 ± 0.08 


1.609 ± 0.061 


2.384 ± 0.177 


NGC 3783 


6.15 ± 0.14 


1.665 ± 0.028 


1.768 ± 0.200 


NGC 4151 


10.11 ± 0.19 


1.519 ± 0.037 


3.466 ± 0.264 


NGC 5548 


4.85 ± 0.13 


1.723 ± 0.028 


0.724 ± 0.149 


NGC 7469 


3.24 ± 0.05 


1.845 ± 0.034 


0.642 ± 0.096 


3C 390.3 


3.44 ± 0.08 


1.644 ± 0.044 


0.453 ± 0.066 


MCG-6-30-15 I 


5.10 ± 0.14 


2.002 ± 0.080 


1.788 ± 0.253 


MCG-6-30-15 II 


5.41 ± 0.10 


2.014 ± 0.038 


1.801 ± 0.226 


NGC 3516 


4.31 ± 0.08 


1.583 ± 0.061 


1.953 ± 0.269 


NGC 3783 


7.56 ± 0.12 


1.689 ± 0.037 


2.080 ± 0.202 


NGC 4151 


22.23 ± 0.24 


1.616 ± 0.055 


4.878 ± 0.795 


NGC 5548 


4.12 ± 0.10 


1.711 ± 0.041 


0.750 ± 0.104 



Note. — The short time scale data are listed in the top half of the table; the bottom 
half lists the long time scale data. Listed are the mean values of -^2-10! F and Fxa- 
Errors for F and Fpca were derived using the formulation of §3.3. The i^2-io errors 
listed above are the average of the errors for F2-10 in a given time bin, which were 
derived from the 2-10 keV continuum light curve, using the mean flux error on the 
~1 ks exposures in that time bin. 
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Table 6 

Fyar VALUES 



Source 


Time 


Continuum 


Line 


Ncinie 


Scale 


^ var \^^J 


R.„^ (%) 

var V ' ^/ 


3C 390.3 


Short 


29.3 ± 0.6 


Undcf. 




Long 


26.9 ± 0.6 


13.7 ± 4.7 


MCG-6-30-15 


Short I 


18.5 ± 0.1 


9.4 ± 2.5 




Short II 


17.4 ± 0.4 


14.4 ± 4.4 




Long I 


16.1 ± 0.7 


14.1 ± 4.8 




T ^ ^ TT 

Long 11 


7.0 ± 0.5 


Undei. 


NGC 3516 


Short 


24.6 ± 0.4 


5.6 ± 2.2 




Long 


32.9 ± 0.4 


14.9 ± 3.7 


NGC 3783 


Short 


4.7 ± 0.7 


6.0 ± 5.5 




Long 


10.0 ± 0.4 


1.5 ± 9.9 


NGC 4151 


Short 


16.7 ± 0.4 


3.2 ± 3.1 




Long 


30.5 ± 0.2 


24.0 ± 4.3 


NGC 5548 


Short 


24.8 ± 0.8 


8.2 ± 12.6 




Long 


39.0 ± 0.7 


28.8 ± 3.9 


NGC 7469 


Short 


13.3 ± 0.3 


Undef. 



Note. — Undefined fractional variability measurements 
indicate a measured variance that is smaller than that ex- 
pected solely from measurement noise. 



Table 7 

F2-io-Fjra Cross-Correlation Results 



Source 


Time 


DCF 


DCF 


IGF 


IGF 


Name 


Scale 


r(d) 


^max 


r(d) 


^max 


3C 390.3 


Short 


-t-12.0 


0.46 


+9+20.2 


0.61 


MCG-6-30-15 I 




+1.3 


0.51 


+1.4+1.5 


0.53 


MCG-6-30-15 II 




-14.4 


0.38 


-14.4+18.0 


0.38 


NGC 3516 




-1-6.4 


0.17 


+6.4+41.2 


0.18 


NGC 3783 




-1-3.3 


0.59 


+4.1+5.4 


0.72 


NGC 4151 




+2.8 


0.23 


+3.5+10.0 


0.27 


NGC 5548 




+6.4 


0.63 


+8.0+6.1 


0.68 


NGC 7469 







0.57 


+0.6+7.2 


0.67 


3C 390.3 


Long 


+51.9 


0.48 


+25.9+159 


0.48 


MCG-6-30-15 I 







0.75 


0+65.6 


0.73 


MCG-6-30-15 II 




-386 


0.45 


0+606 


0.45 


NGC 3516 




-161 


0.66 


-134+262 


0.71 


NGC 3783 







0.53 


0+388 


0.51 


NGC 4151 







0.56 


0+258 


0.61 


NGC 5548 







0.84 


0+428 


0.87 



Note. — Lags are defined such that a positive lag means the continuum 
light curve leads the line light curve. 
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Table 8 

F2_io-r Cross-Correlation Results 



Source 


Time 


DCF 


DCF 


IGF 


IGF 


Name 


Scale 


r(d) 




r(d) 


^TThttX 


3C 390.3 


Short 





0.91 


-3.0±7.8 


0.92 


MCG 6-30-15 I 







0.82 


0±1.4 


0.82 


MCG 6-30-15 II 







0.54 


-0.8±10.9 


0.57 


NGC 3516 







0.90 


0±10.6 


0.90 


NGC 3783 







0.88 


0±5.7 


0.89 


NGC 4151 







0.69 


0±8.4 


0.73 


NGC 5548 




-9.6 


0.85 


-8.0±5.3 


0.88 


NGC 7469 




-3.8 


0.64 


-4.4±8.4 


0.73 


3C 390.3 


Long 





0.85 


0±136 


0.82 


MCG-6-30-15 I 







0.94 


0±66.6 


0.86 


MCG-6-30-15 II 







0.62 


0±555 


0.64 


NGC 3516 







0.92 


0±161 


0.93 


NGC 3783 







0.81 


0±357 


0.81 


NGC 4151 







0.31 


-h21.8±191 


0.46 


NGC 5548 







0.74 


-52.1±506 


0.75 



Note. — Lags are defined such that a positive lag means the continuum 
light curve leads the F light curve. 



Table 9 

r Fa,, Cross-Correlation Resflts 



Source 


Time 


DCF 


DCF 


ICF 


ICF 


Name 


Scale 


r(d) 


^rnax 


T(d) 


^max 


3C 390.3 


Short 


-hl2.0 


0.64 


+9.0+18.1 


0.66 


MCG-6-30-15 I 




+1.5 


0.49 


+1.4+2.0 


0.54 


MCG 6-30-15 II 




+9.6 


0.38 


+10.4+18.2 


0.44 


NGC 3516 




-32.0 


0.25 


-35.2+42.9 


0.33 


NGC 3783 




+3.3 


0.70 


+4.1+5.6 


0.74 


NGC 4151 




+4.2 


0.28 


+4.9+11.6 


0.33 


NGC 5548 




+9.6 


0.43 


+9.6+11.3 


0.47 


NGC 7469 




+5.7 


0.40 


+5.7+9.6 


0.45 


3C 390.3 


Long 


+51.9 


0.52 


+25.9+156 


0.56 


MCG-6-30-15 I 







0.67 


-11.7+97.0 


0.63 


MGG-6-30-15 II 




-386 


0.47 


-386+667 


0.44 


NGC 3516 




-107 


0.79 


-134+300 


0.82 


NGC 3783 







0.47 


0+390 


0.46 


NGC 4151 




-43.6 


0.19 


-174+280 


0.17 


NGC 5548 







0.55 


+52.1+481 


0.61 



Note. — Lags are defined such that a positive lag means the first light 
curve leads the second. 
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Table 10 

Fyar SPECTRAL SAMPLING PARAMETERS 



Source 


MJD Date 




Time Bin 


No. 


Name 


Range 




Size (d) 


Pts 


3C 120 


50458.5-50515, 


.3 


2.0 


29 


3C 390.3 


50220.6-50278, 


.0 


2.0 


29 


Akn 564 


51694.8-51726, 


.7 


1.6 


20 


Fairall 9 


52144.9-52179, 


.0 


1.6 


21 


IC 4329a 


50665.8-50723, 


.5 


2.0 


28 


MCG-6-30-15 I 


50664.1 50672, 


.5 


0.13 


63 


MCG 6-30-15 II 


51622.7 51676, 


.0 


1.6 


34 


NGC 3516 


50523.0-50657, 


.4 


4.8 


28 


NGC 3783 


51960.1-51980, 


.0 


0.66 


30 


NGC 4151 


51870.6 51904, 


.9 


1.2 


29 


NGC 5548 


52091.6 52126, 


.7 


1.3 


27 


NGC 7469 


50244.0 50275 




1.1 


30 


3C 120 


50812.4 51259, 


.7 


20.4 


19 


3C 390.3 


51259.7-51676, 


.0 


20.8 


20 


Akn 120 


50868.1-51259, 


.7 


16.3 


19 


Akn 564 


51678.7-52683, 


.7 


50.2 


20 


Fairall 9 


51679.8-52678, 


.2 


50.0 


20 


MCG- 6-30-15 I 


51926.7 52252, 


.9 


13.1 


23 


MGG-6-30-15 II 


50188.9 51259, 


.7 


53.6 


20 


NGC 3516 


50523.0 51259, 


.7 


28.3 


26 


NGC 3783 


51679.7-52375, 


.9 


29.0 


24 


NGC 4051 


51676.0-52228, 


.3 


27.6 


20 


NGC 4151 


51259.7-51676, 


.0 


17.3 


24 


NGC 5548 


51678.7-52681, 


.4 


45.6 


22 



18 



Sy 1 X-ray Spectral Variability 



3C 390.3 



MCG-6-30-15 I 



1.7 



0.5 



0.2 



220 240 260 

Julian Date (-2450000.5) 



2 5 



3 
2.2 
2 

H 

1.8 



a 

fin 



2 - 



t 

664 



1 I I 



666 668 670 672 
Julian Date (-2450000.5) 



MCG-6-30-15 II 



NGC 3516 



2 5 



3 
2.2 

2 

1.8 

1.6 



B 



+ ++ 



1640 



1660 



Julian Date (-2450000.5) 



1.8 

1.6 

1.4 
3 



H \ 1 1 h 



H \ 1 1 h 



B 



2 



550 



600 



650 



Julian Date (-2450000.5) 



Fig. 1. — Short time scale light curves for continuum flux, T and Fe Ka line flux. The continuum flux is in units of 10 photons cm ^ 
s~^. The Fe flux is in units of 10~* photons cm~^ s~^. The dotted lines respresent the mean values of each parameter. 
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Fig. 1. — Figure 1, cont'd. 
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Fig. 2. — Long time scale light curves for continuum flux, F and Fe Kc< line flux. The continuum flux is in units of 10 photons cm ^ 
~^ . The Fe flux is in units of 10~* photons cm~^ s~^. The dotted lines respresent the mean values of ea<;h parameter. 
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Fig. 3. — Cross-correlation functions for the short time scale light curves. Lags are defined such that a positive lag means the 
curve leads the second. 
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Fig. 5. — Short time scale Fyar spectra. The downward arrows at 6.4 keV indicate the depression in broadband continuum variability due 
to relatively reduced Fe Kq line variability. 
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Fig. 5. — Figure 5, cont'd 
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Fig. 6. — Long time scale Fvar spectra. The downward arrows at 6.4 keV indicate the depression in broadband continuum variability due 
to relatively reduced Fe line variability. 



Markowitz, Edelson & Vaughan 



29 



MCG-6-30-15 II 



NGC 3516 



18 



> 



12 



:+ 1 • 


1 


: +++^ 






5 10 
Energy (keV) 
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